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ABSTRACT 

If  we  want  to  control  vortex  breakdown  on  a 
delta  wing  of  a  military  aircraft  in  order  to  create 
rolling  moments  or  improve  stability  at  high  angle 
of  attack,  the  use  of  CFD  simulations  may  help  us 
to  improve  our  understanding  of  the  flow 
mechanisms  that  occur  with  the  use  of  specific 
control  devices.  The  first  step  of  such  studies  is  to 
make  sure  that  we  can  simulate  properly  and  with 
accuracy  vortex  breakdown  on  a  delta  wing  without 
any  control  device. 

At  Dassault  Aviation  we  have  done  some 
computations  (both  Euler  and  Navier- Stokes)  on 
Mirage  2000  delta  wing  in  a  high  angle  of  attack 
configuration  (with  deflection  of  the  2  slats).  On 
such  a  complex  industrial  configuration  the  flow 
over  the  wing  is  a  mixing  of  the  vortices  coming 
from  the  apex,  the  inner  slat  and  the  outer  slat. 
Euler  computations  lead  to  unsteady  convergence, 
but  the  average  forces  allowed  a  correct  prediction 
of  the  occurrence  of  the  maximum  of  lift  coefficient 
versus  incidence.  Navier-Stokes  computations  for 
the  maximum  lift  coefficient  showed  that  the  wing 
vortices  were  fully  burst,  though  it  was  not  the  case 
on  the  experimental  pressure  distributions.  Several 
ways  of  improvement  are  considered.  The  first  one 
is  to  greatly  improve  the  mesh  discretization  in  the 
region  of  the  vortex.  A  second  one  is  to  improve  the 
turbulence  model:  one  knows  very  well  that 
classical  k-e  models  are  unable  to  simulate  properly 
turbulence  in  rotation,  which  is  the  case  in  the  shear 
layers  that  build  the  vortex  core. 

In  order  to  simplify  the  flow  configuration,  we 
will  show  CFD  results  and  comparison  with  the 
experimental  data  on  a  generic  delta  wing  of  70° 
sweep  angle.  The  corresponding  experiments  were 
conducted  by  ONERA  at  Fauga  F2  wind  tunnel 
[Ref.  1].  We  will  compare  the  vortex  breakdown 
location,  surface  pressure  distributions  and  velocity 
profiles  in  the  vortex  at  locations  upstream  and 
downstream  the  vortex  breakdown.  The  influence 
of  mesh  refinement  in  the  vortex  core  and  of 


turbulence  model  (using  an  EARSM  model  that 
takes  into  account  the  influence  of  rotation  and 
anisotropy)  will  be  presented.  Finally  we  will  give 
some  recommandations  for  the  simulation  of  such 
complex  flows. 

RESUME 

Pour  controler  l'eclatement  tourbillonnaire  sur 
l'aile  delta  d'un  avion  militaire  dans  le  but  de 
generer  des  moments  de  roulis  ou  d'ameliorer  la 
stabilite  a  grande  incidence,  les  simulations 
numeriques  d'ecoulements  nous  aident  a  ameliorer 
la  comprehension  des  topologies  d'ecoulement 
relatives  a  l'utilisation  de  techniques  de  controle 
specifiques.  Une  premiere  etape  dans  ce  genre 
d'etude  est  de  s'assurer  que  l'on  peut  simuler  deja  de 
fagon  precise  l'eclatement  tourbillonnaire  sur  une 
aile  delta  en  l'absence  de  dispositif  de  controle. 

Dassault  Aviation  a  realise  des  calculs  (Euler  et 
Navier-Stokes)  sur  la  configuration  grande 
incidence  (avec  braquage  des  bees)  de  l'aile  delta 
du  Mirage  2000.  Sur  cette  configuration 
industrielle  complexe,  l'ecoulement  sur  l'aile  resulte 
du  melange  des  tourbillons  provenant  de  l'apex,  du 
bee  interne  et  du  bee  externe.  Les  calculs  Euler 
conduisent  a  une  convergence  instationnaire,  mais 
les  efforts  moyennes  predisent  correctement 
l'incidence  de  portance  maximale.  Les  calculs 
Navier-Stokes  realises  pour  l'incidence  de  portance 
maximale  montrent  que  les  tourbillons  de  voilure 
sont  completement  eclates,  alors  que  ce  n'est  pas  le 
cas  sur  les  repartitions  de  pression  experimentales. 
Plusieurs  voies  d' ameliorations  sont  envisagees:  la 
premiere  est  d'enrichir  le  maillage  dans  la  zone  des 
tourbillons.  La  seconde  est  d'ameliorer  la 
modelisation  turbulente:  il  est  connu  que  les 
modeles  k-e  classiques  modelisent  mal  la 
turbulence  en  rotation,  ce  qui  est  le  cas  dans  les 
couches  de  cisaillement  qui  sont  a  l'origine  du 
noyau  visqueux  du  tourbillon. 

De  sorte  a  simplifier  la  configuration 
d'ecoulement,  une  comparaison  calculs  /  experience 
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est  realisee  sur  une  aile  delta  generique  d'angle  de 
fleche  70°.  Les  essais  correspondants  ont  ete 
realises  par  l'ONERA  a  la  soufflerie  F2  du  Fauga 
[Ref.  1].  On  compare  la  position  de  l'eclatement 
tourbillonnaire,  les  repartitions  de  pression  a  la 
surface  de  l'aile,  ainsi  que  les  profils  de  vitesse  dans 
le  tourbillon  en  amont  et  en  aval  de  l'eclatement 
tourbillonnaire.  F'influence  du  raffinement  de 
maillage  dans  le  cceur  du  tourbillon  et  du  modele  de 
turbulence  (utilisation  du  modele  EARSM  qui 
prend  en  compte  l'influence  de  la  rotation  et 
l'anisotropie)  sont  presentes.  En  conclusion, 
certaines  recommandations  sont  donnees  pour  la 
simulation  des  ecoulements  tourbillonnaires. 

1.  INTRODUCTION 

Le  controle  d'un  avion  militaire  a  grande 
incidence,  dans  le  but  d'augmenter  son  agilite, 
necessite  d'une  part  de  pouvoir  controler  les 
tourbillons  de  pointe  avant,  dans  le  but  de  stabiliser 
l'avion  en  contrant  les  dissymetries  tourbillonnaires 
et  aussi  de  disposer  d'une  nouvelle  gouverne  en 
lacet,  d'autre  part  de  pouvoir  controler  l'eclatement 
des  tourbillons  de  voilure  de  fa9on  a  ameliorer  la 
stabilite  longitudinale  et  laterale.  Le  fait  de 
repousser  les  eclatements  tourbillonnaires  a  des 
incidences  plus  elevees  permet,  par  action  sur  le 
Clp,  de  conserver  un  CnP  dynamique  positif. 
D'autre  part,  on  peut  ainsi  ameliorer  la 
caracteristique  Cm  (Cz)  de  l'avion.  En  outre,  par 
une  action  dissymetrique  sur  l'eclatement  des 
tourbillons  droit  et  gauche,  en  anticipant  l'un  et 
retardant  l'autre,  on  peut  generer  des  moments  de 
roulis  importants.  Differents  modes  d'action  ont  ete 
proposes  dans  la  bibliographie:  certains  dispositifs 
sont  de  nature  mecanique,  comme  les  braquages  de 
bees  proposes  par  Rao  [Ref.  2];  d'autres  de  nature 
pneumatique,  comme  le  soufflage  tangentiel  de 
bord  d'attaque  qui  a  ete  teste  par  Wood  [Ref.  3], 
mais  qui  reste  encore  prohibitif  vis-a-vis  du  debit 
de  soufflage  requis.  De  nombreux  autres  dispositifs 
non  cites  ici  figurent  dans  la  bibliographie  et  font 
actuellement  encore  l'objet  d'etudes. 

L'utilisation  de  simulations  numeriques 
d'ecoulement  (CFD)  est  un  outil  d'analyse 
necessaire  a  la  comprehension  de  faction  des 
dispositifs  sur  les  tourbillons  et  leur  eclatement, 
ainsi  qu'a  l'amelioration  de  leur  conception.  Un  pre- 
requis  consiste  a  s'assurer  que  ces  outils  de 
simulation  permettent  de  restituer  fidelement 
l'eclatement  tourbillonnaire  sur  la  voilure  en 
l'absence  de  tout  dispositif  de  controle. 

L'objectif  de  cet  article  est  de  montrer 
1'evaluation  des  codes  de  Dassault  Aviation  vis-a- 
vis  de  problemes  d'eclatement  tourbillonnaire,  tout 
d'abord  sur  une  configuration  de  voilure  delta 
d'avion  militaire,  la  voilure  du  Mirage  2000,  ou 
l'ecoulement  complexe  d'extrados  est  le  resultat  de 
la  coalescence  de  plusieurs  tourbillons,  puis  sur  une 
configuration  plus  simple  d'aile  delta  a  bord 


d'attaque  biseaute  ou  des  donnees  experimentales 
sont  disponibles,  a  la  fois  en  paroi,  ainsi  que  dans  le 
cceur  du  tourbillon. 

2.  ECLATEMENT  TOURBILLONNAIRE  SUR 
LA  VOILURE  DU  MIRAGE  2000: 

La  voilure  du  Mirage  2000  est  un  exemple 
typique  de  configuration  industrielle  qui  est  le  siege 
d'ecoulements  tourbillonnaires  complexes  a  grande 
incidence,  puisque  l'ecoulement  a  l'extrados  resulte 
du  melange  du  tourbillon  d'apex,  de  bee  interne,  de 
bee  exteme,  ainsi  que  de  saumon.  La  voilure  du 
Mirage  2000  a  une  fleche  d'environ  60°,  et  pour  la 
configuration  grandes  incidences  etudiee,  les  bees 
sont  braques  au  maximum. 

Differents  essais  ont  ete  realises  en  soufflerie 
par  l'ONERA  afin  d'identifier  les  efforts  sur  le 
Mirage  2000  en  configuration  grandes  incidences, 
ainsi  que  les  repartitions  de  pression  sur  la  voilure. 
Les  etudes  specifiques  sur  les  tourbillons  de  voilure 
ont  ete  menees  en  soufflerie  sur  une  configuration 
avec  fuselage  generique  a  avant  corps  emousse,  de 
fagon  a  eliminer  de  possibles  interactions  avec  les 
tourbillons  de  pointe  avant.  Une  vue  de  la  maquette 
est  representee  sur  la  figure  1.  Ce  fuselage 
generique  a  ete  defini  de  fa9on  a  etre  representatif 
au  niveau  de  l'apex  de  la  configuration  avion. 

2.1  Calculs  Euler: 

Des  calculs  Euler  ont  ete  menes  sur  l'aile  seule. 
D'autres  avaient  ete  realises  precedemment  sur 
avion  complet.  Le  code  Euler  utilise  est  le  code 
EUGENIE  de  Dassault  Aviation  [Ref.  4].  II  utilise 
des  maillages  non  structures.  Ces  calculs  a  basse 
vitesse  (Mach  0.2)  ont  ete  realises  a  l'aide  d'un 
schema  centre  de  type  Lax-Wendroff,  de  fa9on  a 
limiter  la  dissipation  numerique.  La  figure  2  donne 
un  aper9u  du  maillage  de  l'aile  utilise.  II  est 
important  de  bien  modeliser  la  fente  entre 
l'emplanture  et  le  bee  interne,  puisqu'elle  participe  a 
la  generation  du  tourbillon  majoritaire  de  bee 
interne.  Six  incidences  (al,  a2,  a3,  a4,a5,  a6)  ont 
ete  calculees  de  part  et  d'autre  de  l'incidence  de 
portance  maximale.  La  figure  3  montre 
l'instationnarite  de  l'effort  vertical  Czl  calcule  pour 
ces  differentes  incidences,  cette  instationnarite  etant 
correlee  a  des  battements  des  tourbillons.  Nous 
avons  reporte  la  valeur  moyennee  de  la  contribution 
de  voilure  sur  le  coefficient  de  portance  Cz  (tiree 
des  calculs  Euler)  en  fonction  de  l'incidence  sur  la 
figure  4,  ou  nous  l'avons  comparee  a  la  polaire 
Cz(a)  de  l'avion  complet  tiree  des  essais  en 
soufflerie.  On  constate  que  l'incidence  de  portance 
maximale  est  comprise  entre  a3  et  a4  et  est  assez 
bien  restituee  par  le  calcul  Euler  sur  la  voilure 
seule.  L'ecart  entre  les  points  calcules  et  la  courbe 
experimental  represente  la  contribution  de 
fuselage  au  Cz,  contribution  qui  avait  ete 
precedemment  evaluee  a  partir  des  calculs  sur  avion 
complet. 
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Pour  l'incidence  a3  proche  de  l'incidence  de 
portance  maximale,  nous  avons  compare  les 
repartitions  de  pression  sur  la  voilure  avec  les 
donnees  experimentales  disponibles.  Les  resultats 
sont  representes  selon  differents  arcs  de  cercles 
correspondant  aux  rayons  adimensionnes  suivants 
r/c:  0.34,  0.52,  0.66,  0.75  et  0.95.  Ils  sont 
representes  sur  l'instrumentation  extrados  de  la 
maquette  (cf.  figure  5).  Les  repartitions  de  Kp  selon 
ces  rayons  font  l'objet  de  la  figure  6.  On  y  constate 
que  la  solution  instantannee  finale  du  calcul  Euler 
surestime  la  depression  generee  par  le  tourbillon 
pres  du  bord  d'attaque  et  cree  une  recompression  un 
peu  trop  forte  pres  de  l'emplanture.  Le  tourbillon 
calcule  est  un  peu  plus  intense  et  moins  developpe 
que  le  tourbillon  experimental.  En  moyenne,  par 
effet  de  compensation  entre  la  region  de  bord 
d'attaque  et  la  region  proche  de  l'emplanture,  il 
parait  plausible  que  le  Cz  calcule  soit  proche  du  Cz 
experimental.  II  aurait  ete  plus  juste  de  comparer 
une  moyenne  des  Kp  calcules  aux  resultats 
experimentaux,  mais  celle-ci  n'etait  pas  disponible. 

2.2  Calcul  Navier-Stokes  k-e  en  lois  de  paroi: 

Les  modelisations  fluide  parfait  Euler  ne 
suffisent  plus  des  lors  que  l'on  cherche  a  modeliser 
faction  de  dispositifs,  par  exemple  fluidiques,  sur 
les  couches  limites  avant  qu'elles  ne  decollent  pour 
generer  les  tourbillons.  Les  calculs  Euler 
modelisent  mal  les  effets  d'entrainement.  Aussi,  il 
est  necessaire  d'avoir  recours  a  une  modelisation 
Navier-Stokes  avec  modele  de  turbulence.  Une 
premiere  etape  consiste  alors  a  evaluer  d'abord  ce 
type  de  modelisation  sur  l'aile  sans  dispositif  de 
controle.  Aussi,  nous  avons  calcule  la  configuration 
precedente  du  Mirage  2000  avec  le  code  de  calcul 
Navier-Stokes  AETHER  de  Dassault  Aviation  en 
utilisant  un  modele  de  turbulence  k-e  en 
formulation  lois  de  paroi  [Ref..  5].  L'approche  de 
calcul  utilisee  est  toujours  non  structuree,  ce  qui 
permet  de  faciliter  la  generation  du  maillage. 
L'utilisation  de  la  formulation  en  lois  de  paroi 
permet  d'augmenter  la  taille  de  premiere  maille 
dans  la  couche  limite,  ce  qui,  a  nombre  de  nceuds 
fixe  pour  le  maillage,  permet  de  discretiser  plus 
finement  la  zone  tourbillonnaire,  par  rapport  a  une 
approche  ou  l'on  modeliserait  plus  finement 
l'ecoulement  de  proche  paroi.  La  taille  de  premiere 
maille  retenue  ici  est  de  100  pm.  Le  maillage  non 
structure  correspondant  comporte  environ  993  500 
nceuds.  La  figure  7  illustre  le  maillage  de  la  surface 
de  l'aile,  ainsi  que  celui  du  plan  de  symetrie,  qui  a 
ete  positionne  a  la  jonction  voilure  /  fuselage.  On  y 
constate  que  le  maillage  a  l'extrados  de  la  voilure 
est  assez  regulier  et  assez  dense  pour  pouvoir  capter 
les  tourbillons  avec  precision. 

Un  calcul  Navier-Stokes  k-e  en  lois  de  paroi  a 
ete  mene  pour  l'incidence  a3,  legerement  inferieure 
a  l'incidence  de  portance  maximale.  La 
comparaison  des  Kp  calcules  avec  les  resultats 
experimentaux  selon  les  5  arcs  de  cercle  de  la 


figure  5  fait  l'objet  de  la  figure  8.  Le  calcul  donne 
pour  les  5  sections  des  profils  d'extrados  plats,  qui 
sont  caracteristiques  d'un  tourbillon  eclate  depuis 
l'apex.  Or,  les  donnees  experimentales  montrent 
bien  que  le  tourbillon  resultant  exerce  une 
depression  assez  importante  sur  l'ensemble  de  l'aile, 
y  compris  jusqu'a  la  demiere  section  r/c=0.95,  ce 
qui  contredit  le  resultat  du  calcul.  Il  est  assez 
difficile  de  preciser  l'occurrence  du  point 
d'eclatement,  lorsqu'il  y  a  coalescence,  comme  pour 
cette  configuration  complexe,  de  multiples 
tourbillons.  Un  critere  d'eclatement  a  partir  des 
distributions  de  pression  reste  a  etablir  et  nous  y 
reviendrons  dans  le  cas  de  l'aile  delta  a  bord 
d'attaque  biseaute.  La  figure  9  montre  failure  des 
tourbillons  a  l'extrados  de  l'aile  par  visualisation  de 
la  pression  totale  tiree  du  calcul  k-e  en  lois  de  paroi. 
La  croissance  brusque  des  coeurs  tourbillonnaires 
des  l'apex  est  le  signe  d'un  eclatement  des  f  origine. 

Pour  remedier  a  cet  eclatement  premature  predit 
par  le  calcul  Navier-Stokes  di verses  pistes  sont 
envisageables: 

-  La  premiere  consiste  a  ameliorer  la  modelisation 
turbulente  pour  remedier  a  la  deficience  bien 
connue  du  k-8  qui  modelise  mal  la  turbulence  en 
rotation,  ainsi  que  les  effets  d'anisotropie  de  la 
turbulence.  Un  premier  element  de  reponse  peut 
etre  obtenu  en  utilisant  un  modele  de  turbulence  un 
plus  elabore,  tel  que  l'EARSM  (Explicit  Algebraic 
Reynolds  Stress  Model)  qui  prend  en  compte  les 
effets  de  rotation  et  d'anisotropie  [Ref  6].  Nous 
verrons  au  paragraphe  suivant  le  calcul 
correspondant. 

-  La  seconde  consiste  a  enrichir  le  maillage  dans  la 
zone  du  tourbillon  jusqu'a  obtenir  une  convergence 
de  la  solution  vis-a-vis  d'enrichissements  successifs 
du  maillage.  Cependant,  cette  approche  s'avere 
actuellement  trop  couteuse  en  nombre  de  nceuds 
(done  en  temps  de  calcul)  sur  une  configuration 
aussi  complexe  que  la  voilure  du  Mirage  2000. 
C'est  pourquoi,  l'influence  du  raffinement  de 
maillage  sera  presente  dans  la  partie  3  sur  une 
configuration  d'aile  delta  plus  simple,  celle  de  l'aile 
de  fleche  70°  a  bord  d'attaque  biseaute  qui  a  fait 
l'objet  d'essais  par  l'ONERA  a  la  soufflerie  F2  du 
Fauga. 

2.3  Calcul  Navier-Stokes  EARSM  en  lois  de 
paroi: 

Le  modele  de  turbulence  EARSM  en 
formulation  lois  de  paroi  que  nous  avons  utilise  est 
celui  de  Walin-Johanson.  Nous  avons  recalcule  la 
configuration  precedente  pour  l'incidence  a3  sur  le 
meme  maillage.  La  comparaison  calcul  /  experience 
de  la  figure  10  concemant  les  Kp  dans  les  5 
sections  circulates  precedentes  montre  que  le 
calcul  EARSM  a  ameliore  la  prise  en  compte  de  la 
depression  generee  par  le  tourbillon  dans  les 
premieres  sections.  Si  l'on  retrouve 
approximativement  un  niveau  comparable  a 
l'experience  pres  de  l'apex,  dans  les  sections 
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suivantes  le  calcul  sous-estime  fortement  la 
depression  du  tourbillon  observee  sur  les  resultats 
experimentaux.  Des  la  section  r/c=0.66  le  calcul  fait 
apparaitre  des  profils  de  pression  relativement  plats, 
caracteristiques  d'un  eclatement.  La  visualisation 
des  tourbillons  du  calcul  EARSM  de  la  figure  11 
montre  que  le  tourbillon  de  bee  interne  est  plus 
intense  que  pour  le  calcul  k-e  en  lois  de  paroi  et 
qu'il  coalesce  plus  tard  avec  le  tourbillon  d'apex. 
L'eclatement  du  tourbillon  resultant  est  assez 
difficile  a  reperer.  Cependant  l'affaiblissement  du 
tourbillon  semble  se  produire  au  milieu  du  bee 
interne.  Une  superposition  du  maillage  et  de 
l'e  volution  de  la  pression  totale  dans  cette  zone 
permet  d'apprecier  la  discretisation  utilisee 
localement  pour  decrire  le  tourbillon  (cf.  figure  12). 
On  y  constate  que  le  tourbillon  est  decrit 
diametralement  par  une  vingtaine  de  mailles,  ce  qui 
s'avere  deja  une  discretisation  raisonnable.  Les 
effets  d'influence  de  maillage  vont  etre  abordes  sur 
la  configuration  suivante. 

3.  ECLATEMENT  TOURBILLONNAIRE  SUR 
L’AILE  DELTA  DE  70°  DE  FLECHE 

L'interet  principal  de  calculer  cette  nouvelle 
configuration,  outre  le  fait  qu'elle  est  plus  simple, 
e'est  qu'elle  est  assez  riche  en  mesures 
experimentales,  a  la  fois  en  ce  qui  concerne  les 
repartitions  de  pression  a  l'extrados  et  aussi  les 
mesures  de  vitesses  dans  le  cceur  du  tourbillon  pour 
plusieurs  sections  en  amont  et  en  aval  de 
l'eclatement  tourbillonnaire.  De  plus,  l'etude 
experimentale  a  ete  menee  dans  une  plage 
d'incidence  assez  large  allant  jusqu'a  40°.  Nous 
nous  limiterons  cependant  aux  validations  pour  les 
incidences  26°  et  27°,  ou  les  effets  de  confinement 
dans  la  veine  de  F2,  non  represents  dans  les 
calculs,  n'ont  pas  encore  un  grand  effet  sur  la 
position  de  l'eclatement.  Les  conditions  d'essai 
correspondent  a  une  pression  generatrice  Pi  0  de  1 
bar  et  une  temperature  generatrice  de  288  K.  La 
vitesse  soufflerie  est  de  24  m/s  ce  qui  conduit  a  un 
nombre  de  Reynolds  par  metre  de  1,65  million. 

La  configuration  geometrique  est  celle  d'une 
aile  delta  de  70°  de  fleche  a  bord  d'attaque  biseaute. 
Cette  maquette  est  constitute  d'une  plaque  plane  de 
20  mm  d'epaisseur  avec  des  bords  d'attaque 
biseautes  a  15°  cote  intrados.  Sa  corde  d'emplanture 
fait  950  mm.  Une  visualisation  de  la  geometrie  de 
la  demi-voilure  est  montree  sur  la  figure  13 
representant  les  divers  maillages  de  surface  utilises. 

3.1  Maillages  utilises  pour  cette  etude: 

Un  maillage  grossier  preliminaire  de  cette 
configuration  a  ete  realise  dans  un  premier  temps 
afin  de  pouvoir  identifier  la  zone  ou  se  developpe  le 
tourbillon  et  pouvoir  reenrichir  celle-ci  sur  un 
second  maillage  discretise  beaucoup  plus  finement. 
La  figure  13  montre  le  maillage  peau  de  la  demi- 
voilure  pour  le  maillage  preliminaire  et  le  maillage 
raffine.  Sur  ce  dernier,  le  biseau  a  ete  raffine  de 


fagon  a  disposer  de  beaucoup  de  precision  dans  la 
zone  du  contournement  du  bord  d'attaque.  De  plus, 
on  a  fait  en  sorte  d'avoir  des  tailles  de  mailles  de 
l'ordre  du  millimetre  dans  toute  la  zone  de  la 
voilure  ou  s'exerce  la  depression  du  tourbillon.  En 
ce  qui  concerne  le  maillage  espace,  une 
discretisation  des  premieres  mailles  compatible  d'un 
calcul  en  lois  de  paroi  a  ete  appliquee.  Par  contre, 
les  tailles  de  mailles  obtenues  dans  l'espace  sont 
pilotees  par  les  tailles  des  mailles  a  la  peau.  La 
figure  14  montre  la  comparaison  des  maillages  dans 
le  plan  de  coupe  X=700  mm  et  illustre  les 
differences  de  densites  de  mailles  evoquees 
precedemment.  Le  maillage  preliminaire  comporte 
143  200  nceuds,  alors  que  le  maillage  raffine 
comporte  1  119  300  nceuds  environ. 

L'influence  de  la  discretisation  se  voit 
directement  sur  la  simulation  de  l'eclatement 
tourbillonnaire.  La  figure  15  montre  pour  ces  deux 
maillages  la  visualisation  du  cceur  tourbillonnaire 
au  dessus  de  l'aile  delta  par  l'intermediaire  d'une  iso 
valeur  de  la  pression  totale.  On  constate  que  le 
calcul  sur  maillage  grossier  presente  un  tourbillon 
mou  qui  n'eclate  pas,  alors  que  le  maillage  fin  fait 
apparaitre  un  tourbillon  plus  intense  dont  le  noyau 
grossit  brusquement  a  l'endroit  du  point 
d'eclatement. 

3.2  Caracterisation  du  tourbillon: 

Sur  le  maillage  fin,  deux  calculs  ont  ete  menes 
pour  les  incidences  de  26°  (ou  l'on  dispose  des  Kp 
experimentaux)  et  de  27°  (ou  l'on  dispose  des 
mesures  de  vitesse  dans  le  tourbillon).  Deux 
modeles  de  turbulence  ont  ete  testes:  le  k-8  en  lois 
de  paroi  et  le  modele  EARSM  en  lois  de  paroi.  Sur 
la  configuration  a  26°,  nous  avons  repere  la  position 
de  l'axe  du  tourbillon  issu  du  calcul  EARSM.  On 
obtient  une  fleche  de  l'axe  de  75°  et  une  incidence 
de  l'axe  de  7°  par  rapport  a  la  voilure  (cf.  figure  16). 
Nous  avons  compare  ces  valeurs  aux  correlations 
empiriques  etablies  par  Erickson  [Ref.  7],  qui 
donnent  une  fleche  de  l'axe  tourbillonnaire  de  77° 
et  une  incidence  de  l'axe  de  8°.  On  obtient  des 
valeurs  assez  proches. 

En  ce  qui  concerne  la  position  de  l'eclatement, 
elle  a  ete  deduite  des  calculs  en  determinant 
l'endroit  ou  la  vitesse  axiale  sur  l'axe  du  tourbillon 
s'annule.  Les  positions  de  l'eclatement 
tourbillonnaire  tirees  de  ces  calculs  ont  ete 
reportees  sur  la  figure  17  qui  montre  pour  les 
differentes  incidences  etudiees  en  essais  les 
positions  extremales  reperees  pour  l'eclatement.  Les 
differences  observees  en  essais  sur  les  releves  de  la 
position  d'eclatement  sont  partiellement  lies  au  fait 
que  certaines  mesures  correspondent  a  des 
balayages  croissant  ou  decroissant  en  incidence, 
d'autre  part  des  differences  ont  ete  observees  sur  les 
points  d'eclatement  entre  cote  droit  et  cote  gauche 
de  la  maquette.  Pour  les  deux  incidences  calculees, 
on  constate  que  le  calcul  avec  modele  k-s  predit  un 
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eclatement  un  peu  plus  precoce  que  le  calcul  avec 
modele  EARSM.  Cependant,  les  deux  calculs 
donnent  une  position  assez  proche  de  la  position  la 
plus  aval  tiree  des  essais. 

3.3  Repartition  de  pression  a  l’extrados  de  la 
maquette: 

Les  comparaisons  avec  l'experience  sont 
effectuees  pour  les  calculs  a  l'incidence  de  26°  . 
Nous  avons  etudie  l'influence  du  maillage  et  celle 
du  modele  de  turbulence. 

En  ce  qui  concerne  l'influence  de  maillage,  la 
figure  18  montre  les  distributions  de  Kp  a  l'extrados 
pour  les  deux  maillages  decrits  precedemment.  Sur 
le  maillage  preliminaire  les  tourbillons  assez  mous 
generent  une  depression  peu  intense.  Sur  le 
maillage  raffine,  les  depressions  sont  bien  plus 
intenses  et  l'on  y  observe  deux  caracteristiques  de 
l'eclatement  tourbillonnaire:  L'eclatement  etant 
situe  pour  ce  calcul  EARSM  a  l'abscisse  700  mm, 
d'une  part  la  zone  de  depression  est  pincee  a 
l'endroit  de  l'eclatement,  puis  s'elargit,  d'autre  part 
le  tourbillon  devient  plus  mou  en  aval  de 
l'eclatement  se  qui  se  traduit  par  une  diminution  du 
gradient  de  pression  axial  apres  eclatement  (les 
courbes  de  Kp  selon  l'envergure  evoluent  moins 
apres  eclatement).  La  figure  19  compare  a 
l'experience  les  niveaux  de  Kp  obtenus  avec  les 
deux  maillages  pour  les  sections  X=400  mm  (avant 
eclatement)  et  X=800  mm  (apres  eclatement).  Le 
maillage  fin  permet  de  restituer  les  pics  de 
depression,  contrairement  au  maillage  grossier, 
dont  les  tourbillons  ne  sont  pas  assez  intenses. 

L'effet  du  modele  de  turbulence  employe  est 
represente  sur  la  figure  20  pour  les  deux  memes 
sections  dans  le  cas  du  maillage  fin.  On  note  assez 
peu  de  difference  sur  la  prediction  des  niveaux  de 
depression,  mis  a  part  que  la  zone  de  depression  est 
un  peu  plus  large  pour  le  k-8  apres  eclatement,  etant 
donne  qu'il  fait  eclater  plus  tot  le  tourbillon  que 
l'EARSM. 

La  comparaison  avec  les  mesures 
experimentales  montre  un  assez  bon  accord, 
notamment  avec  l'EARSM,  apres  eclatement  du 
tourbillon  pour  le  plan  X=800  mm.  En  ce  qui 
concerne  la  section  avant  eclatement  a  X=400  mm, 
des  differences  apparaissent  sur  la  largeur  des  pics 
de  depression  et  sur  les  niveaux  de  pression  pres  du 
bord  d'attaque.  En  effet,  les  calculs  sont  realises  en 
turbulent,  alors  qu'experimentalement  a  cette 
abscisse  sur  la  maquette  l'ecoulement  est 
probablement  encore  laminaire,  ce  qui  a  pour  effet 
de  faire  apparaitre  un  tourbillon  secondaire 
correspondant  au  plateau  de  Kp  pres  du  bord 
d'attaque.  Ce  tourbillon  secondaire  a  aussi  pour 
effet  de  decaler  le  tourbillon  primaire  legerement 
vers  l'emplanture. 


3.4  Validation  des  vitesses  au  coeur  du 
tourbillon: 

Les  comparaisons  avec  l'experience  sont 
effectuees  pour  les  calculs  a  l'incidence  de  27°  . 
Nous  avons  etudie  l'influence  du  modele  de 
turbulence.  Nous  ne  presenterons  que  les  vitesses 
pour  des  sondages  en  envergure  passant  par  l'axe  du 
tourbillons,  les  sondages  obtenus  a  d'autres 
positions  verticales  presentant  des  caracteristiques 
similaires.  Les  comparaisons  sont  presentees,  etant 
donne  les  mesures  disponibles,  en  vitesse 
longitudinale  Vx  et  verticale  Vz,  et  non  pas  en 
composantes  axiale  et  orthoradiale.  Cependant 
l'ecart  entre  ces  deux  representations  est  faible. 

Sur  la  figure  21,  on  montre  la  comparaison  des 
Vx  avant  et  apres  eclatement.  En  amont  de 
l'eclatement,  on  note  peu  de  difference  entre  les 
deux  modeles  de  turbulence,  sauf  a  partir  de 
l'abscisse  X=600  mm,  etant  donne  que  le  k-8 
conduit  a  un  eclatement  plus  precoce  que  le  modele 
EARSM.  La  largeur  du  pic  de  Vx  est  bien  predite, 
mais  le  maximum  de  vitesse  longitudinale  est  sous- 
estime.  Apres  eclatement,  la  bulle  de  recirculation 
est  moins  large  pour  l'experience  que  pour  les 
calculs  et  son  profil  de  type  sillage  se  comble  plus 
vite.  Le  calcul  k-8  presente  un  sillage  qui  se  remplit 
plus  vite  que  le  calcul  EARSM,  etant  donne  qu'il 
presente  un  eclatement  plus  precoce. 

La  figure  22  illustre  la  composante  Vz  avant  et 
apres  eclatement.  On  note  peu  de  differences  entre 
les  deux  modeles  de  turbulence.  Cependant,  les 
deux  calculs  font  apparaitre  un  cisaillement  moins 
fort  au  cceur  du  noyau  tourbillonnaire  que  les 
resultats  experimentaux.  Deux  explications 
possibles  a  cela:  soit  la  discretisation  du  noyau 
tourbillonnaire  est  insuffisante,  soit  le  modele  de 
turbulence  conduit  a  un  noyau  tourbillonnaire  trop 
visqueux.  Apres  eclatement,  le  cisaillement  sur  la 
composante  Vz  diminue  avec  l'intensite  du 
tourbillon. 

3.5  Examen  du  maillage  et  recommandations: 

L'utilisation  d'un  maillage  tres  raffine  au  niveau 
de  la  peau  a  permis  de  generer  un  maillage  a 
mailles  isotropes  de  dimension  caracteristique  de 
l'ordre  de  1  mm  a  proximite  de  la  surface  de  l'aile. 
Cependant,  au  dela  d'une  hauteur  de  50  mm,  pour 
limiter  le  nombre  de  nceuds,  la  hauteur  des  mailles 
augmente,  et,  du  fait  des  condensations,  la  largeur 
des  mailles  augmente  aussi.  Les  tailles  de  mailles 
deviennent  alors  de  l'ordre  de  3  a  4  mm,  notamment 
dans  la  region  ou  se  produit  l'eclatement 
tourbillonnaire,  comme  le  fait  apparaitre  la  figure 
23  qui  presente  le  maillage  dans  le  plan  vertical 
passant  par  l'axe  du  tourbillon.  Un  tel  maillage  est 
relativement  fin  et  precis  pour  capter  l'origine  du 
tourbillon.  Ce  tourbillon  se  convecte  ensuite  et 
eclate  dans  une  zone  ou  le  maillage  est  moins 
dense. 
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Ce  qu'il  faudrait  etudier,  c'est  la  convergence  de 
la  solution  vis-a-vis  d'enrichissements  successifs 
dans  le  cceur  du  tourbillon,  afin  de  determiner  des 
criteres  portant  sur  le  nombre  de  points  necessaires 
a  repartir  selon  un  diametre  pour  une  discretisation 
precise  du  cceur  tourbillonnaire.  Un  tel  probleme  de 
sous-maillage  pour  un  calcul  Navier-Stokes  s'avere 
non  trivial  si  l'on  desire  ne  pas  augmenter  de 
maniere  prohibitive  le  nombre  de  nceuds  du 
maillage  espace.  Les  criteres  de  detection  du  cceur 
tourbillonnaire  peuvent  etre  deduits  de  la  pression 
totale  ou  bien  de  la  composante  axiale  du 
rotationnel.  Une  partie  du  probleme  consiste  a 
isoler  la  zone  de  proche  paroi,  afin  de  ne  pas  y 
propager  les  raffinements. 

4.  CONCLUSION 

La  simulation  numerique  de  l'eclatement 
tourbillonnaire  sur  aile  delta  a  ete  etudiee  par 
Dassault  Aviation  sur  deux  configurations 
distinctes: 

-  la  configuration  industrielle  de  la  voilure  du 
Mirage  2000  avec  braquage  de  bees  pour  les 
grandes  incidences,  ou  l'ecoulement  d'extrados 
resulte  de  la  coalescence  de  plusieurs  tourbillons, 

-  la  configuration  fondamentale  de  l'aile  delta  de 
70°  de  fleche  qui  a  ete  bien  identifiee  en  soufflerie 
vis-a-vis  de  l'eclatement  tourbillonnaire. 

La  modelisation  Euler  donne  des  niveaux  de 
portance  assez  representatifs.  En  revanche,  les 
calculs  Navier-Stokes  sous-estiment  celle-ci  pour  la 
voilure  du  Mirage  2000.  L'approche  EARSM 
ameliore  un  peu  la  modelisation  de  l'intensite  des 
tourbillons  par  rapport  a  un  calcul  k-8. 

L'utilisation  d'un  maillage  tres  fin  a  permis  de 
representer  de  fagon  tres  satisfaisante  les 
repartitions  de  pression  sur  l'aile  delta  de  l'ONERA. 
Cependant  le  pic  de  vitesse  axiale  et  le  cisaillement 
de  vitesse  orthoradiale  sont  sous-estimes  par  les 
modeles  de  turbulence  encore  trop  visqueux. 

Des  efforts  restent  a  faire  pour  ameliorer  la  prise  en 
compte  des  effets  de  rotation  dans  les  modeles 
turbulents,  ainsi  que  pour  pouvoir  enrichir  de  fagon 
astucieuse  et  a  un  cout  non  prohibitif  en  nombre  de 
nceuds  les  zones  de  noyaux  tourbillonnaires  sur  des 
configurations  industrielles  d'ailes  delta  d'avions 
militaires. 
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Figure  1:  maquette  de  la  voilure  du  Mirage  2000 
avec  fuselage  generique  (maquette  ONERA) 


Maillage  de  1'exlrados  de  la  voilure  du  Mirage  2000  pour  le  calcul  EUGENIE 
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Figure  2:  maillage  de  la  surface  de  l'aile  du  Mirage 
2000  pour  les  calculs  EUGENIE 
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Convergence  de  i'effort  vertical  Czl  pour  les  6  incidences 

Figure  3:  convergence  en  effort  normal  des  calculs 
Euler 


Figure  4:  Comparaison  en  Cz  entre  calculs  Euler 
sur  la  voilure  et  Cz  experimental  sur  avion  complet 


Figure  5:  Instrumentation  de  la  maquette  de  voilure 
du  Mirage  2000 


Figure  7 :  maillage  Navier-Stokes  de  l'aile  du 
Mirage  2000 


Figure  8:  repartition  des  Kp  selon  les  5  sections  en 
envergure  pour  le  calcul  k-8  a  l'incidence  a3 


Aile  Mirage  2000  bees  braques  -  a  3  -  calcul  Navier-Stokes  K-eps  lois  de  paroi 


Figure  9:  tourbillons  d'extrados  du  calcul  k-8 


Figure  6:  repartition  des  Kp  selon  les  5  sections  en 
envergure  pour  le  calcul  Euler  a  l'incidence  a3 


Figure  10:  repartition  des  Kp  selon  les  5  sections  en 
envergure  pour  le  calcul  EARSM  a  l'incidence  a3 
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Aile  Mirage  2000  bees  braques  -  a3  -  calcul  Navier-Stokes  EARSM  lois  de  paroi 


Figure  11:  tourbillons  d'extrados  du  calcul  EARSM 


Figure  12:  Discretisation  du  tourbillon  de  bee 
interne  avant  eclatement 


Figure  13:  Visualisation  des  maillages  peau  de  la 
demi-aile  delta  pour  le  maillage  preliminaire  et  le 
maillage  raffine 


Figure  14:  Coupes  dans  le  plan  X=700  mm  pour  les 
deux  maillages 


Figure  15:  Visualisation  du  tourbillon  par  la 
pression  totale  sur  le  maillage  preliminaire,  puis  sur 
le  maillage  raffine 


0  100  200  300  400  500  600  700  800  900  1000  1100 

X 


Figure  16:  Position  de  l'axe  tourbillonnaire  en  vue 
de  dessus  pour  les  calculs  k-8  et  EARSM 

Aile  Delta  avec  bord  d'attaque  aigu  J  —  £  £SE  “t? 

O  calcul  k-e  lol  de  parol 

Mach=0.071  (V=24m/s)  «  calcul  EARSM  lol  de  parol 


Figure  17:  Comparaison  calculs  /  essais  sur  la 
position  d'eclatement  du  tourbillon 
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Coupe  X=800mm 
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Figure  19:  Comparaison  des  Kp  pour  les  deux 
maillages  avant  et  apres  eclatement 


Figure  20:  Comparaison  des  Kp  pour  les  modeles 
k-8  et  EARSM  avant  et  apres  eclatement 


Figure  21:  Comparaison  de  1'evolution  du  Vx  en 
envergure  dans  le  cceur  du  tourbillon  entre  calculs 
et  essais 


Figure  22:  Comparaison  de  1'evolution  du  Vz  en 
envergure  dans  le  cceur  du  tourbillon  entre  calculs 
et  essais 
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Figure  23:  Discretisation  du  tourbillon  pres  du  point 
d'eclatement 
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Paper:  23 

Author:  Mr.  Rosenblum 

Question  by  Mr.  Sacher:  Excellent  piece  of  work.  You  choose  the  most  challenging  case, 
the  delta  wing  with  a  deflected  leading-edge  flap.  Question:  Did  you  first  compute  the  case 
with  plain  configuration  and  observe  better  agreement  between  CFD  vs.  EFD? 

Answer:  Yes,  we  computed  the  case  with  the  plain  configuration  and  observed  the  same  sort 
of  differences  between  k-e  and  EARSM  turbulence  models.  But  we  didn’t  have  experimental 
data  to  compare  with,  since  the  slot/deflected  configuration  is  the  one  for  high  angle  of  attack. 

Question  by  Dr.  Hummel:  What  is  the  reason  for  the  differences  in  the  vortex  breakdown 
location  in  the  F2-wind  tunnel?  Are  the  results  not  reproducible? 

Answer:  As  it  was  mentioned  by  Mr.  Molton  in  Paper  19  and  by  Dr.  Mitchell  in  Paper  20,  the 
two  points  on  the  experimental  curve  correspond  to  the  limits  of  the  fluctuation  of  the  location 
of  vortex  breakdown. 


